PROCEEDINGS 

°F SCIENCE 



oo 

o 
o 

> 
O 

(N 



iu: 

£^ 

(N 
> 

in 

00 

o 
o 

d 

oo 

o 



Topological susceptibility in (2+1)-flavor lattice QCD 
with overlap fermion 



T.W. Chiu* 'JS. Aoki^S S. Hashimoto^^, T.H. Hsieh«, T. Kaneko^^, H. Matsufuru*^ 
J. NoakP, T. Onogi'', N. Yamada^^ (for JLQCD and TWQCD Collaborations) 

" Physics Department, Center for Theoretical Sciences, and Center for Quantum Science and 
Engineering, National Taiwan University, Taipei 10617, Taiwan 

^ Graduate School of Pure and Applied Sciences, University of Tsukuba, Tsukuba 305-8571, 
Japan 

Riken BNL Research Center, Brookhaven National Laboratory, Upton, NY11973, USA 
" High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801, Japan 
f School of High Energy Accelerator Science, The Graduate University for Advanced Studies 

( Sokendai), Tsukuba 305-0801, Japan 
^ Research Center for Applied Sciences, Academia Sinica, Taipei 115, Taiwan 
^ Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 



We determine the topological susceptibility Xt in the topologically-trivial sector generated by 
lattice simulations of Nf = 2 + I QCD with overlap Dirac fermion, on a 16^ x 48 lattice with 
lattice spacing 0.11 fm, for five sea quark masses m^/ ranging from OTj /6 to m^ (where m^ is 
the physical strange quark mass). The Xt is extracted from the plateau (at large time separation) 
of the 2-point and 4-point time-correlation functions of the flavor-singlet pseudoscalar meson Tj', 
which arises from the finite size effect due to fixed topology. In the small nig regime, our result 



of Xt agrees with the chiral effective theory. Using the formula Xt = ^('Wh +"^^/ 
Leutwyler-Smilga, we obtain the chiral condensate I.^^{2 GeV) = [249(4)(2)MeV]^ 



by 



The XXVI International Symposium on Lattice Field Theory 
July 14-19 2008 
Williamsburg, Virginia, USA 



* Speaker. 

^ Email: twchiu@phys.ntu.edu.tw 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Topological susceptibility in (2+1 )-flavor lattice QCD 



T.W. Chiu 



1. Introduction 

In Quantum Chromodynamics (QCD), the topological susceptibility (Xt) is the most crucial 
quantity to measure the topological charge fluctuations of the QCD vacuum, which plays an im- 
portant role in breaking the Ua{1) symmetry. Theoretically, Xt is defined as 

Xt= J d\{p{x)p{0)) , p(x) = ^^£^,vACTtr[F^vWf)icTW] (1.1) 

where p (x) is the topological charge density expressed in term of the matrix-valued field tensor 
F^y. With mild assumptions, Witten [|l[ and Veneziano obtained a relationship between the 
topological susceptibility in the quenched approximation and the mass of rj' meson (flavor singlet) 
in full QCD with Nf degenerate flavors, namely, (quenched) = flm^,/{4Nf) where /^r = 131 
MeV, the decay constant of pion. This implies that the mass of T]' is essentially due to the axial 
anomaly relating to non-trivial topological charge fluctuations, which can turn out to be nonzero 
even in the chiral limit, unlike those of the (non-singlet) approximate Goldstone bosons. 

Using the Chiral Perturbation Theory (ChPT), Leutwyler and Smilga obtained the follow- 
ing relation in the chiral limit 

Xt = ^^ ^ -X + ^i'^l)^ iNf = 2 + \), (1.2) 

— + — + — 

where m„, m^, and are the quark masses, and £ is the chiral condensate. This implies that in the 
chiral limit (m„ 0) the topological susceptibility is suppressed due to internal quark loops. Most 



importantly, ( |1.2| ) provides a viable way to extract £ from Xt in the chiral limit. 



From (1.1), one obtains 



Xt = ^, Qt^jd'xp{x), 

where Q. is the volume of the system, and Qt is the topological charge. Thus, one can obtain Xt by 
counting the number of gauge configurations for each topological sector. Obviously, for a set of 
gauge configurations with Qt = 0, it gives Xt = 0- However, even for a topologically-trivial gauge 
configuration, it may possess non-trivial topological excitations in sub-volumes. Thus, one can 
measure Xt using the correlation of the topological charges of two sub-volumes. 

In general, for any topological sector with Qt, using saddle point expansion on the QCD par- 
tition function in a finite volume, it can be shown that [Q] 

lim {p{x)pm = l^(^-Xt-^^+m-'), (1-3) 



V // Q.\Q. IXt^ 

where C4 = — ^ [(Gf )e=o — '^{Qi)'^e=o\ ■ However, for lattice QCD, it is difficult to extract p(x) and 
Qt unambiguously from the gauge link variables, due to their rather strong fluctuations. 
To circumvent this difficulty, one may consider the Atiyah-Singer index theorem ^ 

Qt = n+-n^=mdex{&), (1.4) 

where n± is the number of zero modes of the massless Dirac operator '3 = y^ {d^ + igA^ ) with ± 
chirality. Since ^ is anti-Hermitian and chirally symmetric, its nonzero eigenmodes must come 
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in complex conjugate pairs (i.e., ^0 = /A0 implies ^75^ = —/A 750, for A = A* 7^ 0) with zero 
chirality (/ d^x^^Y5^ = 0). Thus one can obtain the identity 

„,-„_=/</^,™,rte(S. + ».)-'(x,,)l, (1.5) 

by spectral decomposition, where the nonzero modes drop out due to zero chirality. In view of ( p^ 
and (1.5), one can regard m ti{Y5{& + m)^^ {x,x)] as topological charge density, to replace p{x) in 
the measurement of Xt- 

For lattice QCD, it is well-known that the overlap Dirac operator ||6|, 0] in a topologically 
non-trivial gauge background possesses exact zero modes (with definite chirality) satisfying the 
Atiyah-Singer index theorem. Writing the massive overlap Dirac operator as 

D(m)= Uo + -) + Uo--)75- °^ 



2^ ' ir^^jHi^y 

where mo) is the standard Hermitian Wilson operator with negative mass —mo (0 < mo < 2), 
then the topological charge density can be defined as 

p^(x) =mtr[75(D, + m)-]], 

where (D^. + m)^^ is the valence quark propagator with quark mass m, and Dc is a chirally sym- 
metric operator relating to D(0) by Dc = D{0)[1 - D(0)/(2mo)]"' Here p^ix) is justified to 
be topological charge density, since it can be shown that Y,xPm{x) = ?i+ — which is similar to 
its counterpart in continuum, ([L^). 

Now we can replace p{x) with Pi(x), and use (13) to extract Xt for any topological sector. 
However, on a finite lattice, it is contaminated by mji, m,,/ and any states which can couple to 
(pi (jc)pi (0)). An alternative is to consider the correlator of the flavor-singlet pseudoscalar meson 

17' ii 

which is equal to the disconnected part {p\^{x])p\{x2)) q at large separation, but it tends to the 
asymptotic value faster than the later since it only couples to the states containing r]' . Then the 
time-correlation function of Tj' is fitted to A + B{e^'^' + e^'^^'^^'^) to obtain the constant A = 

\ \ (Q] _ 



j' ^"^^ ^^"^ which to extract Xt provided that |c4| ^ ^Xt^- This was how 
we determined the topological susceptibility in two-flavor lattice QCD with fixed topology p^ ]. 

However, it was unclear to what extent the assumption |c4| <C ^xfO. was satisfied. To elim- 
inate this constraint, we compute the 4-point correlator of f]' , as well as the 2-point correlator. 
Theoretically, in a fixed topology, the former behaves as |Q] 

^^^Wm)- ■ ^ I'Mh = f (1 - ^ + + ^(^-"''"'-"') + m-'), (1.7) 

From ( L6) and (1/7), one can solve for Xt and C4 (or equivalently, the parameter y) 



Xt = ^ + ^[2k2-Vh/3), (1.8) 



, = ^ = -i^^^l^fl-^) (19) 
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where —k2 and ^4 are the asymptotic values of 2-point and 4-point correlators at large separation. 
It is interesting to note that if one neglects the y term in (1^) and ([LTI), they reduce to 



(1.10) 

(1.11) 



which provide two independent estimates of Xt- In other words, if |y| ^ 1, then ( [Ol ), ( |1.10D and 
( |1.11[ ) all give compatible results for Xt- On the other hand, if ( |1.10| ) and ( 1.11 ) turn out to be 
quite different from each other, then \y \ must be substantially larger than zero, and a more reliable 
estimate of Xt could be given by ( |1.1 1[ ). 

For the (2+l)-flavor QCD, the rj' interpolating operator must take into account of the fact that 
different flavors have different quark masses, namely. 



Nf 



Nf 



m, 



ff=i 



Nf 



(1.12) 



where rjj is called the "topological" rj' operator for computing topological charge density correla- 
tors. 

In this paper, we use 80 pairs of low-lying eigenmodes of the overlap operator to evaluate the 
2-point and 4-point correlators of Tjf , and to extract their asympototic values —^2 and ^4. Then we 



use (1.8)-(1.S) to obtain Xt and y. Note that C4 is related to the leading anomalous contribution to 
the rj' — rj' scattering amplitude in QCD, as well as the dependence of the vacuum energy on the 
vacuum angle 6. 
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Figure 1: Low-mode saturation of (a) the 2-point function C^/_ (f ) (b) the 4-point function C4^^ (?) 



2. Lattice Setup 



Our simulations are carried out in the topologically-trivial sector (with Qt = 0) for (2-1-1)- 
flavor QCD on a 16^ x 48 lattice at a lattice spacing 0. 1 1 fm (for details, see [11] and [p^]). For 
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the gluon part, the Iwasaki action is used at j8 = 2.30, together with unphysical Wilson fermions 



and associated twisted-mass ghosts jl3|]. The unphysical degrees of freedom generate a factor 
det[//^.(— mo)/(//^(— mo) + /i^)] in the partition function (we take mo = 1.6 and jj. = 0.2) that 
suppresses the near-zero eigenvalue of //„,(— mo) and thus makes the numerical operation with the 
overlap operator substantially faster. Furthermore, since the exact zero eigenvalue is forbidden, the 
global topological change is preserved during the molecular dynamics evolution of the gauge field. 

For m^ = 0.100, we take five sea quark mass m„(^) values: 0.015, 0.025, 0.035, 0.050, and 
0.100 that cover the mass range m^/6-mj. After discarding 500 trajectories for thermalization, 
we accumulate 2500 trajectories in total for each sea quark mass. In the calculation of Xt, we 
take one configuration every 5 trajectories, thus we have 500 configurations for each niq. For 
each configuration, 80 pairs of lowest-lying eigenmodes of the overlap-Dirac operator D(0) are 
calculated using the implicitly restarted Lanczos algorithm and stored for the later use. 

3. Results 

In practice, we use 80 pairs of low-lying eigenmodes of the overlap operator to evaluate the 
2-point and 4-point time-correlation functions of rjj 

1 ^ 1 
Cr^^it) = 73^LL(T]r(^2,i< + 0^f(^b"))> Jim73C^f(0 = -^2, 

C4^/_(0 = -j-^^(T]r(x4,M + 3f)T7y(x3,M + 2f)T7y(f2,M + 0^f(^b"))> Jim TgQr^f (0 = ^4- 

»=1 .Y; 

Thus it is crucial to check whether these 80 eigenmodes suffice to saturate Cfj/,(f) and C4f,/,(?) re- 
spectively. In Fig. H, we plot Cf,/, (t) and C^j^i^ (t) for m„ = 0.015, versus the number of eigenmodes 
(nev) 20, 40, 60, and 80 respectively. Obviously, C^/, (t) is well saturated with 80 eigenmodes for 
the time range 15 < f < 24 where it attains a plateau. Similarly, C4^!^{t) is also well saturated for 
the time range 9 < ? < 14 where it attains a plateau. The low-mode saturation also holds for all five 
sea quark masses. 

In Fig. ^ we plot the values of Xt'^'^ (1-8) and y ([T^) versus the sea quark mass m^a, together 



with the values of Xt obtained from the 2-point function ( |1.10| ) and the 4-point function ( |l . 1 1| ) 
respectively. Evidently, the values of Xt from (L8), ( 1.10 ), and ( |1.11[ ) are in good agreement with 



one another. For the smallest four quark masses, 0.015, 0.025, 0.035, and 0.050, the data points of 
a'^Xt are well fitted by the ChPT formula [|] 

Xt = , ^, _i , (3.1) 

m„ + m^ + Ms 

with a^r = 0.0021(1). In order to convert £ to that in the MS scheme, we calculate the renor- 
malization factor (2 GeV) using the non-perturbative renormalization technique through the 
RI/MOM scheme. Our result is Z^{{mathrm2 GeV) = 0.833(8) 0. With a-^ = 1833(12) 
MeV determined with ro = 0.49 fm [|l2|], the value of £ is transcribed to 

1^(2 GeV) = [249(4) (2)MeV]^ 
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Figure 2: Topological susceptibility Xr^'^ and y = C4/{2x}^) versus sea quark mass niqa for (2+l)-flavor 
lattice QCD with fixed topological charge Qt — 0. 
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Figure 3: The topological susceptibility Xt versus niq for lattice QCD with fixed topology Qt — 0. 

which is in good agreement that extracted from Xt = {Q})/^ with Qt determined by the spectral 
flow method for the 2+1 flavors QCD configurations generated by the RBC and UKQCD Collab- 
orations with domain-wall fermions [p3|]. Also, it is in good agreement with our previous results 
extracted from Xt in 2-flavor QCD [T^], and in the £-regime from the low-lying eigenvalues 
[p^. The errors represent a combined statistical error (a^^ and Z^^) and the systematic error 
respectively. 

At this point, it is instructive to plot Xt versus m^, for 2-flavor QCD (data from H^), and 
(2-i-l)-flavor QCD (this work), as shown in Fig. |3[ Now we can see clearly how the topological 
susceptibility changes with respect to the number of flavors. 

4. Concluding remark 

In this paper, we have obtained the topological susceptibility Xt and C4 in (2-1-1 )-flavor QCD 
from a lattice calculation of 2-point and 4-point correlators at a fixed global topological charge 
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Qt = 0. The expected sea quark mass dependence of Xt from ChPT is clearly observed. Our result 
asserts that the topologically non-trivial excitations are in fact locally active in the QCD vacuum, 
even when the global topological charge is zero. The information of these topological excitations 
is carried by the low-lying eigenmodes of the overlap Dirac operator. We will use the values of 
Xt we have determined to remove the artifacts due to the fixed topology in a finite volume and to 



obtain the physical results in the 6 vacuum 1 17 



Finally we note that our result of the ratio \c4\/Xt is substantially less than one in the chiral 
limit, similar to its counterpart in quenched QCD [jTsj, [T^. This seems to suggest that the quantum 
corrections would suppress the emergence of dilute instanton gas in the full QCD vacuum. 

Numerical simulations are performed on Hitachi SRI 1000 and IBM System Blue Gene Solu- 
tion at High Energy Accelerator Research Organization (KEK) under a support of its Large Scale 
Simulation Program (No. 08-05), and also on IBM and HP clusters at NCHC and NTU-CC in Tai- 
wan. This work is supported in part by the Grant-in-Aid of the Japanese Ministry of Education 
(Nos. 18340075, 18740167, 19540286, 19740160, 20025010, 20039005, 20340047, 20740156), 
the National Science Council of Taiwan (Nos. NSC96-2112-M-002-020-MY3, NSC96-2112-M- 
001-017-MY3, NSC97-2119-M-002-001), andNTU-CQSE (Nos. 97R0066-65, 97R0066-69). 
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